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THE DARK MATTER RADIAL PROFILE IN THE CORE OF THE RELAXED CLUSTER A2589 
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ABSTRACT 

We present an analysis of a Chandra- ACIS observation of the galaxy cluster A2589 to constrain the 
radial distribution of the total gravitating matter and the dark matter in the core of the cluster. A2589 is 
especially well-suited for this analysis because the hot gas in its core region (r < O.lrvir) is undisturbed 
by interactions with a central radio source. From the largest radius probed (r = O.OTr-vir) down to 
r « 0.02rvir dark matter dominates the gravitating mass. Over this region the radial profiles of the 
gravitating and dark matter are fitted well by the NEW and Hernquist profiles predicted by CDM. The 
density profiles are also described well by power laws, p oc r~", where a = 1.37±0.I4 for the gravitating 
matter and a = 1.35 ± 0.21 for the dark matter. These values are consistent with profiles of CDM halos 
but are significantly larger than a « 0.5 found in LSB galaxies and expected from self-interacting dark 
matter models. 

Subject headings: galaxies:clusters:individual (A2589) — dark matter — intergalactic medium — 
X-ray:galaxies:clusters — cosmological parameters 



1. INTRODUCTION 

The very precise constraints that have been placed re- 
cently on the cosmological world model by observations 
of the cosmic microwave background (e.g., Spergel et al. 
2003) and high-redshift supernovae (e.g., Perlmutter et al. 
1999) require that most of the matter in the universe is 
non-baryonic "dark matter" and that the largest contrib- 
utor to the energy density of the universe is the "dark en- 
ergy" . But these two mysterious quantities that dominate 
the energy density of the universe still could be merely fit- 
ting parameters - akin to Ptolemaic epicycles - that are 
not physical quantities. There is, therefore, great urgency 
to discover the nature of the dark matter and dark en- 
ergy which are the foundations of the new cosmological 
paradigm. 

The structure of dark matter (DM) halos is a sensitive 
probe of the properties of the DM. In the standard ACDM 
paradigm, N-body simulations show that the radial den- 
sity profiles of DM halos are fairly uniform, approximately 
parameterized by the NEW profile, p(r) oc r^^{rs + r)^^ 
(Navarro et al. 1997). At small radii (r <C rg) the 
density profiles follow a power law, although the precise 
value of the power-law exponent remains controversial; i.e., 
p{r) oc r~", with a = 1 according to NEW and a = 1.5 
according to the simulations of Moore et al. (1999). Ob- 
servations of the rotation curves of low-surface brightness 
(LSB) galaxies (e.g., Swaters et al. 2000) suggest a pro- 
file for the DM that is substantially flatter {a « 0.5) than 
ACDM in the central regions. These observations inspired 
Spergel & Steinhardt (2000) to propose the existence of 
"self-interacting" DM (SIDM). In the SIDM model the 
DM particles are assumed to possess some cross section 
for elastic collisions with each other. Detailed CDM simu- 
lations incorporating the SIDM idea have confirmed that 
the DM profiles of LSB galaxies can be flattened as ob- 
served (e.g., Dave et al. 2001). 

Like LSB galaxies, galaxy clusters provide excellent 
venues to study DM because they are DM-dominated from 



deep down into their cores (« O.Olrvir; e.g., Dubinski 1998) 
out to their virial radii (Irvir). Moreover, several powerful 
techniques are available for probing the cluster matter dis- 
tributions: optical studies of galaxy dynamics and gravita- 
tional lensing, and X-ray observations of the hot gas. Each 
of these techniques has advantages and disadvantages. For 
example, some advantages of using X-ray observations are 
that the hot gas in clusters traces the three-dimensional 
gravitational potential with an isotropic pressure tensor. 
The quality of data for the hot gas is limited only by the 
sensitivity and resolution of the X-ray detectors and not 
by the flnite number of galaxies. 

The most important limitation associated with X-ray 
observations of cluster mass distributions is the assump- 
tion of hydrostatic equilibrium. A large fraction of nearby 
clusters (z < 0.2) do have regular image morphologies and 
appear to be relaxed on 0.5-1 Mpc scales (e.g., Mohr et al. 
1995; Buote & Tsai 1996; Jones & Forman 1999). For such 
clusters with regular morphologies (i.e., not currently ex- 
periencing a major merger) cosmological N-body simula- 
tions have determined that masses calculated by assuming 
perfect hydrostatic equilibrium are generally quite accu- 
rate (e.g., Tsai et al. 1994; Evrard et al. 1996; Mathiesen 
et al. 1999). But clusters that are observed to be relaxed 
on 0.5-1 Mpc scales are typically associated with cooling 
flows (e.g., Buote & Tsai 1996), and Chandra observa- 
tions have demonstrated that the inner cores of cooling 
flows are highly disturbed exhibiting holes and fllamen- 
tary structures (e.g., Fabian et al. 2000; David et al. 2001; 
Ettori et al. 2002) which certainly raise serious questions 
about the assumption of hydrostatic equilibrium. 

Hence, for studies of DM it is imperative to find clusters 
that are relaxed (i.e., have smooth, regular X-ray images) 
from deep within their cores out to ~ Mpc scales. Our 
recent analysis of the Chandra ACIS-S data of one such 
cluster (A2029, Lewis et al. 2002, 2003) has provided the 
strongest constraint on the core DM profile for a galaxy 
cluster to date; i.e., a — 1.19 ± 0.04, constrained down to 
« O.Olrvir. The core DM density profile of A2029 is con- 
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sistent with the NFW profile and rules out an important 
contribution from SIDM in this cluster. Since the DM 
properties of A2029 may not be typical, and the expected 
cosmological scatter (Bullock et al. 2001) needs to be ad- 
dressed, more X-ray observations of clusters with undis- 
turbed cores are very much needed. 

In our present study we analyze of the core DM profile 
of the galaxy cluster A2589 using a new Chandra ACIS-S 
observation. We selected A2589 as the brightest, nearby 
(z = 0.0414) cluster with a smooth X-ray morphology 
(according to ROSAT images, Buote & Tsai 1996) and 
without a known central radio source (Bauer et al. 2000). 
The latter criterion greatly increases the likelihood that 
the X-ray image of the core is undisturbed. (The redshift 
of A2589 corresponds to an angular diameter distance of 
171 Mpc and 1" = 0.83 kpc assuming Hq = 70 km s~^ 
Mpc-S Om = 0.3, and Oa = 0.7.) 

The paper is organized as follows. In §2 we present the 
observation and discuss the data reduction. The spectral 
analysis of the ACIS-S data is discussed in §3. We cal- 
culate the gravitating matter distribution in §4, the gas 
mass and gas fraction in §5, and estimate the stellar mass 
in §6. The DM profile and the systematic error budget 
are discussed in §7 and §8 respectively. Finally, in §9 we 
present our conclusions. 

2. OBSERVATION AND DATA REDUCTION 

A2589 was observed with the ACIS-S CCD camera for 
approximately 15 ks during AO-3 as part of the Chan- 
dra Guest Observer Program. The events list was cor- 
rected for charge-transfer inefficiency according to Towns- 
ley et al. (2002), and only events characterized by the stan- 
dard ASCA grades^ were used. The exposure time after 
applying these procedures is 13.7 ks. The standard CIAO^ 
software (version 2.3) was used for most of the subsequent 
data preparation. 

Since the diffuse X-ray emission of A2589 fills the en- 
tire S3 chip, we attempted to use the standard background 
templates^ to model the background. Although there were 
no strong flares (i.e., "strong" being defined as count rates 
> 10 times nominal), there was mild flaring activity during 
most of the observation. Consequently, after running the 
standard lc_CLEAN script to clean the source events list of 
flares with the same screening criteria as the background 
templates, only 3 ks remain. 

To salvage a large portion of the observation (8.7 ks) 
while obtaining a more accurate estimate for the back- 
ground than given by the standard templates alone, we 
followed the procedure discussed in Buote et al. (2003a). 
That is, we flrst extracted the source spectrum from re- 
gions near the edge of the S3 chip and subtracted from 
it the spectrum of the corresponding regions of the back- 
ground template. (Following Markevitch et al. 2003 we 
renormalized the background template by comparing the 
ACIS-S 1 data of the source and background observations 
for energies above 10 keV. This resulted in a normaliza- 
tion 8% above nominal.) We fitted the resultant spectrum 
with a two-component model consisting of a thermal com- 
ponent, represented by an apec thermal plasma, and a 

^ littp:/ /cxc. harvard.edu/udocs/docs/docs. html 
^ http;//cxc.harvard.edu/ciao/ 
^ http://cxc.haxvaxd.edu/cal 



broken power-law (BPL) model, representing the residual 
flaring background which is most pronounced at high en- 
ergies. For the BPL model we obtain a break energy of 
5 keV with power-law indices of below the break and 1.9 
above the break. We scale this model to the area of the 
annular regions of interest (see below) and generate a pha 
correction file for use in XSPEC; i.e., for each region an- 
alyzed below we subtract background contributions from 
the standard templates and a correction pha flle account- 
ing for the excess flaring. 

3. SPECTRAL ANALYSIS 

In Figure 1 we display the ACIS-S3 image of A2589. 
The isophotes are regularly shaped and are approximately 
elliptical with e = 0.2 — 0.3 over the region of interest. The 
centroid of the innermost contour is displaced with respect 
to that of the outer contours shown by « 15". Very im- 
portantly, we see no evidence for disturbances or holes in 
the hot gas that are ubiquitous in the cores of cooling flow 
clusters as discussed in §1. In this paper we ignore the 
flattening of the X-ray isophotes and instead obtain cir- 
cularly averaged spectral quantities (i.e., gas density and 
temperature) on the sky to construct the spherically aver- 
aged deprojected mass distribution. 

We extracted spectra in concentric circular annuli lo- 
cated at the X-ray centroid computed within a radius of 
20" with initial center on the peak of the X-ray emission. 
(We examine the sensitivity of our results to the chosen 
center position in §8.) The widths of the annuli were cho- 
sen to have approximately equal background-subtracted 
counts in the 0.5-5 keV band. Because of the high back- 
ground we did not use the data for R>3', and the widths 
of the annuli had to be relatively large. These restrictions 
resulted in a set of five annuli within R = 3.2' as listed in 
Table 1. We constructed ARF files for each annulus using 
the latest corrections that account for the time-dependent 
degradation in the quantum efficiency at lower energies. A 
single RMF file appropriate for the CTI-corrected S3 data 
was used. 

For each annulus we fitted the background-subtracted 
spectrum with an APEC thermal plasma modified by 
Galactic absorption (4.15 x 10^° cm^^). For the APEC 
model we usually let the iron abundance be a free pa- 
rameter, and the abundances of all the other elements are 
tied to iron in their solar ratios; i.e., we fit a metallicity. 
The spectral fitting was performed with XSPEC (Arnaud 
1996) using the method. Hence, we rcbinned all spec- 
tral channels to have a minimum of 30 counts per energy 
bin. We take the solar abundances in xspec (vll.2.0bc) 
to be those given by Grevesse & Sauval (1998) which use 
the correct "new" photospheric value for iron which agrees 
also with the value obtained from solar-system meteorites 
(e.g., Mc William 1997). We restricted the spectral fit- 
ting between 0.5-5 keV to avoid calibration uncertainties 
at lower energies and systematic errors associated with 
the background subtraction at high energies. To estimate 
the statistical errors on the fitted parameters we simulated 
spectra for each annulus using the best-fitting models and 
fit the simulated spectra in exactly the same manner as 
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Fig. 1. — Chandra ACIS-S image (0.5-5 keV) of a 4.5' x 4.5' region centered on A2589. The image has been smoothed with a gaussian 
filter with a = 12". Smoothed contours overlaid axe logarithmically spaced in intensity. Note the image has not been corrected for exposure 
variations, nor has the background been subtracted. 



Table 1 

Parameters from the Spectral Fits 









Rout 




T 




norm 




Annulus 


(aJTcmin) 


(kpc) 


(arcmin) 


(kpc) 


(kcV) 


(solar) 


(10-'-^ cm-'') 


(xVdof) 


1 








0.53 


26 


3.16 ± 0.24 


1.19 ±0.35 


1.22 ±0.10 


53.6/56 


2 


0.53 


26 


0.82 


40 


3.14 ±0.24 


0.65 ±0.24 


1.41 ±0.10 


68.7/57 


3 


0.82 


40 


1.34 


66 


3.30 ± 0.20 


0.75 ±0.21 


2.72 ±0.15 


118.5/104 


4 


1.34 


66 


1.85 


91 


3.03 ± 0.19 


0.64 ±0.16 


2.86 ±0.14 


128.2/107 


5 


1.85 


91 


3.18 


156 


3.18 ±0.16 


0.36 ±0.09 


6.74 ±0.21 


165.1/167 



Note. — Results of fitting a single apec plasma emission model modified by Galactic absorption directly to the annular spectra (i.e., without 
spectral deprojection). The norm parameter is the emission measure of the apec model as defined in xspec: 10~^*(/ rieripdV) / in D'^ {1 + 2^) 
with units cm~^. The quoted errors are la computed using the Monte Carlo procedure described in §3. 




done for the actual data. From 100 Monte Carlo simu- 
lations we compute the standard deviation for each free 
parameter which we quote as the "la" error. (We note 

that these Icr error estimates generally agree very well 
with those obtained using the standard Ax^ approach in 

XSPEC.) 

The parameters obtained from the spectral fits are listed 
in Table 1. The simple one-temperature model is a good 
fit to the data in all annuli. In Figure 2 we show the ACIS- 



S3 spectra of the innermost and outermost annuli and the 
associated best-fitting one-temperature model. Visual in- 
spection of Figure 2 reveals that the spectral shapes of 
annuli #1 and ^5 are similar above ^ 2 keV indicating 
that they have similar temperatures. Near 1 keV the spec- 
trum of annulus #1 is more peaked reflecting its larger iron 
abundance. In sum, the one-temperature model reveals a 
nearly isothermal gas but with evidence for a metallicity 
gradient similar to those observed in other clusters with 
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a dominant central galaxy (e.g., Lewis et al. 2002; Ettori 
et al. 2002; Gastaldello & Molendi 2002; Sanders & Fabian 

2002) . 

These results agree with those obtained using ROSAT 
PSPC data by David ct al. (1996) within the rela- 
tively large statistical errors associated with the PSPC 
data. David et al. (1996) obtain a temperature of T = 
1.7to 3 keV (90% conf.) within R = V which is marginally 
inconsistent with our results obtained with Chandra. 
Howciver, spectral parameters obtained from X-ray spec- 
tra of groups and clusters can depend on a variety factors 
such as bandwidth (e.g., see Buote et al. 2003a,b). David 
et al. (1996) fitted models over 0.3-2 kcV appropriate for 
the PSPC. If we use 0.5-2 keV we obtain T « 2.5 kcV 
within R = 30" consistent with the PSPC result. 

Although the fits are acceptable, we investigated 
whether they could be improved further. We allowed other 
elemental abundances to be free parameters. It was found 
that allowing silicon and oxygen to be free did offer sig- 
nificant (though not substantial) improvements in the fits 
in some annuli. Typically, both abundances fitted to low 
values (i.e., sub-solar ratios with respect to iron) but are 
consistent with the iron abundance within the estimated 
Icr errors. Since it was also found that when allowing sil- 
icon and oxygen to be free the fitted temperatures and 
normalizations did not vary appreciably within their la 
errors, we decided to keep silicon and oxygen tied to iron 
in the fits as done for the other elements. 

We also found that adding another temperature com- 
ponent did not improve the fits much. For example, in 
the central annulus adding another temperature compo- 
nent improved by 3 but with the addition of two free 
parameters (one for temperature and one for normaliza- 
tion - abundances are tied to the first temperature com- 
ponent). The second temperature component (T w 1 keV) 
has only « 1% of the emission measure indicating that 
a single temperature component clearly dominates the 
ACIS-S spectra. If we add a high-temperature compo- 
nent (T > 10 keV) to annulus #4 the fit is improved by 
Ax^ = 20. This weak high-temperature component is 
almost certainly associated with incomplete background 
subtraction in this annulus. In §8 we discuss errors asso- 
ciated with the background and other systematic effects. 

4. THE RADIAL PROFILE OF GRAVITATING MASS 

4.1. Method 

The approach we use to calculate the mass distribution 

follows closely our previous study of A2029 (Lewis et al. 

2003) . We approximate the hot gas in the cluster as spher- 
ic-ally symmetric and in hydrostatic equilibrium so that the 
total gravitating mass is, 

kBT{r) ( dlnpg dlnT 



grav(< r) - T ^^^^ 



dhir dhir 



(1) 



where r is the radius in a three-dimensional volume, pg is 
the gas mass density, T is the gas temperature, nip is the 
atomic mass unit, fi is the mean atomic weight of the gas 
(taken to be 0.62), ks is Boltzmann's constant, and G is 
Newton's constant. Mgrav is the enclosed gravitating mass 
which is the sum of the masses of the stars, the gas, and 
the DM. We evaluate the derivatives of pg and T using 
parameterized models as discussed below. We note that 



spherical averaging of the data is appropriate for compari- 
son to the spherically averaged DM profiles obtained from 
cosmological simulations which is a key goal of this paper. 

We note that the cooling time of the gas is « 1 Gyr 
within the central annulus, indicating that some heat 
source is required to prevent the gas from cooling. Mod- 
els of cluster cooling flows including both feedback from 
an AGN and thermal conduction appear to successfully 
prevent the gas from cooling (e.g., Brighenti & Mathews 
2003). In the context of this paradigm, since no substan- 
tial morphological disturbances in the X-ray image are cur- 
rently observed in the core of A2589, this cluster presum- 
ably has had sufficient time to relax since the previous 
feedback episode. 



4.2. Deprojected Radial Profiles of Gas Density and 
Temperature 

We first attempted to obtain three-dimensional gas pa- 
rameters by performing a spectral deprojection procedure 
based on the "onion-peeling" method that we have used 
in previous studies (Buote 2000; Buote et al. 2003a; Lewis 
et al. 2003). Because the S/N of the A2589 data is rela- 
tively low, we arc unable to obtain precise constraints on 
the deprojected temperature and abundance profiles using 
this procedure (but see §8). Consequently, we projected 
parameterized models of the three-dimensional quantities, 
pg and T, and fitted these projected models to the re- 
sults obtained from our analysis of the data projected on 
the sky (Table 1). In this manner we obtained good con- 
straints on the three-dimensional radial profiles of pg and 
T. Note that for each annulus we assign a single radius 
value, r = [{rll^l + rf/^)/2]2/3, where nn and rout are re- 
spectively the inner and outer radii of the annulus /shell. 
We found this prescription to provide an excellent esti- 
mate of the mean weighted radius in our previous study 
of A2029 (Lewis et al. 2003). 

The norm parameter of the apec model obtained from 
the spectral fits (see Table 1) is proportional to / pg^dV, 
where dV represents the emitting volume. For our data in 
projection we have dV — Adl, where A is the area of the 
annulus on the sky and dl is the length element along the 
line of sight; i.e., J Pg^dl oc norm/A. (It is assumed that 
the plasma emissivity is constant along the line of sight 
which is a very good approximation for the nearly isother- 
mal spectrum of A2589.) In Figure 3 we plot / Pg^dl along 
with the best-fitting /? model (Cavaliere & Fusco-Femiano 
1978), 



r(i)r(3/3-i) 



-3/3/2 



Pg dl = 



r(3/3) (l+i?2/r2)3/3-o.5' 



(2) 

(3) 



where the integral is evaluated along the line of sight 
(— oo, -|-oo), and R is the radius of the annulus on the sky.. 
The best-fitting parameters, la errors, and values are 
listed in Table 2. The errors are calculated by fitting the 
norm profiles obtained for each Monte Carlo error simula- 
tion (§3). We find that the simple /? model is an excellent 
fit. 

We modeled the temperature data as the emission- 
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Fig. 3. — Left Panel: Chandra radial profile of the projected gas density squared (J pg^dl) obtained by dividing the norm parameter of the 
APEC model (Table 1) by the area of the annulus as described in §4.2. Horizontal bars indicate the sizes of the annuli used to extract spectra, 
and the sizes of the spherical shells used in our deprojected analysis. Overlaid is the best-fitting /3-model {solid curve). Right Panel: Chandra 
radial temperature profile of A2589. Overlaid are the best-fitting isothermal [dashed (red) curve), and power-law {solid (blue) curve) models. 
Also shown are the la limits for the power-law model {dotted (green) curves). The temperature data points refer to projected quantities 
(Table 1) while the temperature models are deprojected (Table 2). 



Table 2 

Deprojected Radial Profiles of Gas Density and Temperature 





Gas Density 


Model 


(xVdof) 


(kpc) f3 


Peo 

(10^26 g cm-3) 


/3 


0.8/2 


40 ± 12 0.39 ±0.04 


1.93 ± 0.12 



Temperature 







Tioo 




Model 


(xVdof) 


(keV) 


P 


isothermal 


1.0/4 


3.16 ± 0.22 





power law 


1.0/3 


3.17 ±0.43 


-0.006 ± 0.06 



Note. — We find a central electron density of rieg = 9.8 ± 0.6 X 10 cm (the conversion factor between n^Q and pgQ is 5.09 X 10^''). Tioo 
is the temperature of the isothermal model (p = 0), and is the temperature evaluated at r = 100 kpc for the power law model. 



weighted projection of a power law profile, 

p 



T{r) - T] 



100 



100 kpc 



(4) 



(r(i?))- / p^^T, 



dl/ / pg^dl, (5) 



+ P 

f 00 kpc 

where Tioo is the temperature at r = f 00 kpc, pg is the (3 
model given above, and again R is the radius of the annu- 
lus on the sky. The best-fitting model {{T{R))) is plotted 
in Figure 3 and the parameters and la errors for T(r) are 
given in Table 2. The temperature profile is nearly isother- 



mal. This is illustrated by the excellent fit obtained for p 
fixed to (Figure 3 and Table 2). 

4.3. Three- Dimensional Radial Mass Profile 

Using the parameterized functions for the gas density 
and temperature just described, the radial mass distribu- 
tion (eq. 1) becomes. 



Mgrav(< r) = r 



fcsTioo(r/100kpc)P / 3(3ryr, 



Gprrip 



1 + r'^ /r"^ 



P 



(6) 
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Table 3 

Parameters of Gravitating Matter and Dark Matter Radial Profiles 



Mass Model 


(x'/dof) 


Mum "Ts 
7 (1013 Mq) (kpc) 


c (Mpc) 




Gravitating Matter 


power law 
NFW 


5.0/3 
2.0/3 


1.63 ±0.14 1.11 ±0.44 

153 ± 362 


6.6 ±2.5 1.01 ±0.32 



Dark Matter 



power law 2.7/2 1.65 ± 0.21 0.90 ± 0.33 

NFW 1.9/2 ■_ ■_ 208 ±865 4.9 ± 2.4 1.03 ±0.51 



Note. — The power-law temperature model is used for all the fits. The fits to the DM profile exclude the inner mass point and assume 

Afstars/Lv = 6.1 (§6). 




Fig. 4. — Left Panel: Total gravitating mass of the cluster for each radial bin obtained from the isothermal (open squares, red) and the 
power-law (open circles, blue) fits to the gas temperatures. The /3-model parameterization for pg was used in both cases. The mass points 
are fitted with corresponding power law functions {solid line, blue: power-law T model, dashed line, red: isothermal T model). Right Panel: 
Total gravitating mass (data points enclosed with open circles), overlaid with two different fitted models: NFW ((blue) solid curve) and 
power-law ((red) dashed line). Each model uses the power-law parameterization of the temperature profile. The enclosed gas mass is plotted 
as data points enclosed with open triangles. Also plotted is the mass range expected from the stars ((green) dot- dashed curves): The lower 
curve assumes Msturs / Ly = IMq/Lq and the upper curve Mstars/iy = 12Mq/Lq (see §6). We have used large open symbols to identify 
the data points, as some of the error bars are difficult to see in this logarithmic plot. The upper axis shows the radius in units of the virial 
radius, r^ii = 1.72 Mpc. 



In Figure 4 we plot Afg^av evaluated for each radial bin 
obtained from the isothermal (open squares, red) and the 
power-law {open circles, blue) fits to the gas temperatures. 
The mass profiles obtained from the isothermal and power- 
law models agree extremely well. Fitting a power law, 

to the mass profile itself gives 7 — 1.63 ± 0.14 and 
7 = 1.61 ± 0.13 for the power-law (p) and isothermal 



{p — 0) temperature parameterizations respectively. The 
excellent agreement demonstrates that the mass profile is 
not particularly sensitive to the temperature parameteri- 
zation, and thus we focus on the power-law temperature 
parameterization henceforth. 

The NFW density and enclosed mass profile are given 
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by, 



p(r) 



Pc{z)Sc 



(r/rs)(l +r/r,)2 



M(< r) = 47rpc(-2)4r-f 



In 



rs +r 



rs+r 



(8) 



(9) 



where Pc{z) = 3H{z)'^ /SttG is the eritical density of the 
universe at redshift z, Vs is the scale radius, and Sc is 
a characteristic dimensionless density which, when ex- 
pressed in terms of the concentration parameter c, takes 
the form, 

200 



3 ln(l + c) - c/(l + c) ■ 
Hence, there are two free parameters for the NFW mass 
profile: rs and c. The virial radius, defined to be the radius 
where the enclosed average mass density equals 2Q0pc{z), 
is simply ryir = CTg. We calculate Pc{z) at the redshift of 
A2589 assuming = 0.3 and I^a = 0.7 today. 

In Figure 4 we show the result of fitting an NFW model 
to the mass profile along with the power-law fit for com- 
parison. The parameters for both models are listed in 
Table 3. The NFW model provides a very good fit to 
mass points, better than the power law model, though the 
power law is still a fair representation of the mass profile. 
Both models lie above the innermost data point. Exclud- 
ing the innermost data point does not significantly change 
the results. For example, when omitting the inner mass 
point the power-law mass fit gives 7 = 1.47 ± 0.15 which 
is within la of the value obtained when fitting all the 
points (Table 3). Similarly, we obtain a concentration pa- 
rameter c = 7.8 ± 3.4 for NFW when excluding the inner 
point which is quite consistent with the value given in Ta- 
ble 3. (We mention that a Hernquist model (Hernquist 
1990) fits as well as NFW (x^ = 1.9 for 3 dof), but the 
Moore model is not quite as good a fit (x^ = 6.7 for 3 

The virial radius (rvir) obtained from the NFW fit to 
A^grav is rather small (1.01±0.32 Mpc) for a massive galaxy 
cluster, though our value does have a large uncertainty. It 
is consistent with the virial radius, r^ir = 1.72 ±0.32 Mpc, 
obtained using the virial mass temperature relation of 
Mathiesen & Evrard (2001) with the temperature of the 
isothermal model (Table 1). 

5. GAS MASS, GAS FRACTION, AND flm 

The radial profile of the enclosed mass of hot gas ob- 
tained from the /3-model (eq. 2) is plotted in Figure 4. 
The gas fraction (/gas = M^^/Mg,^^) is 0.059 ± 0.033 in 
the inner radial bin and rises to 0.087 ± 0.006 in the last 
radial bin; i.e., the gas mass is a minor contributor to the 
gravitating mass over the region studied. The gas frac- 
tion in the last radial bin may be safely considered to be 
a lower limit since simulations and other measurements at 
larger radii in clusters suggest /gas = {Q.2 — Q.2>)h^Q^^ (e.g., 
Allen et al. 2002). Assuming the global baryon fraction in 
A2589 is representative, then we place an upper limit on 
the present matter density, < f^b//gas = 0.56 ± 0.05 
where we have used n^h"^ = 0.024± 0.001 for h = 0.7 from 
recent measurements of the CMB (Spergel et al. 2003). 
We emphasize that our upper limit on results from un- 
derestimating the baryon fraction because (1) we do not 
measure a global value of /gas, and (2) we account for only 
the baryons associated with the hot gas. 



6. GRAVITATING MASS-TO-STELLAR LIGHT RATIO AND 
THE STELLAR MASS CONTRIBUTION 

To calculate the ratio of gravitating mass to stellar light 
we used the fits to the F-band surface brightness profile of 
A2589 published by Malumuth & Kirshner (1985). These 
authors found that a King model with core radius Tc = 
0.98/7,7Q^ kpc and a total luminosity of 2.3 x IQ^^Lv.qH^q 
provides a a good fit to surface brightness profile out to 
R = I'iQh^Q kpc. The King model fit translates to a vol- 
ume light density profile, /Oy(r) = py(0)(l + r'^/r^)^'^/^, 
where p\/(0) = 3.8-L\/.Qpc^'^. The enclosed luminosity as 
a function of radius in three dimensions {Lv{< r)) is cal- 
culated by integrating pv over the volume of the sphere of 
radius r. We plot Mgi.av(< f)/Lv{< r) in Figure 5. 

-Wgrav/-^y climbs from 2.8±1.5Mq/L0 at the inner ra- 
dial bin to 73.6 ± 4.8 Mq/L© at the outer radial bin. At 
the resolution of our chosen radial bin sizes there is no ev- 
idence for a sharp transition in Mgrav/i'y over the region 
studied. 

Converting Ly to stellar mass (Mgtars) is problem- 
atic because of the uncertain stellar mass-to-light ratio 
{Mstaxs/ Lv) of elliptical galaxies. Often stellar dynam- 
ical measurements of the cores of elliptical galaxies are 
taken as estimates for the stellar mass. But since such 
dynamical studies (like X-rays) are sensitive to the total 
gravitating mass, they cannot be relied upon to yield ro- 
bust values of the stellar mass separated from the DM. 
Stellar population synthesis studies offer a direct method 
to calculate Msturs/ Lv but unfortunately allow for a large 
range of values (Mstars/.^^y = 1 — I2M0/L0) because of 
uncertainties associated with the stellar initial mass func- 
tion, the age of the population(s), and the metallicity (e.g.. 
Pickles 1985; Maraston 1999) ~ though see Mathews'(1989) 
for arguments that the stellar mass is similar to the value 
determined from stellar dynamics. 

In Figure 4 we plot Mstars for the cases Mgtars /i-y = 
1Mq/Lq and Mgtars/ Lv = I2M0/L0 to show the al- 
lowed range from population synthesis studies of ellipti- 
cal galaxies. The upper limit for Mgtars essentially passes 
through the value of Mgrav at the second radial bin (r = 
33 kpc) but greatly exceeds Mgj-av in the central radial bin. 
If Afgrav in the central bin is reliable then Mgtars /-^v must 
be small enough to match MgravZ-^v = 2.8 ± 1.5 Mq/ Lq 
obtained from the X-ray analysis. Hence, we assume a 
plausible range of stellar mass-to- light ratios, Mgtars/ Ly = 
(1 — 7 A)Mq/ Lq^ where the upper limit is the 3(t upper 
limit on Afgrav/^y in the central radial bin and the lower 
limit is determined by the population synthesis models. 



7. THE DARK MATTER RADIAL PROFILE 

Inspection of Figure 4 reveals that the gravitating mass 
exceeds the combined masses of stars and hot gas for 
r > 0.02rvir indicating that DM dominates over most 
of the region studied. We write the mass in DM as, 

MdM = Mgrav - Mgag - Mgtars, whcrC Mgas and Mstars 

are described in §5 and §6. For Afgtars we define a fidu- 
cial value by requiring that interior to the effective radius 
{Re = 32.6 kpc, Malumuth & Kirshner 1985) we have 
Mdm = Mgtars. This definition allows for a fairly smooth 
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Fig. 5. — Ratio of gravitating mass to stellar light from the cD galaxy in A2589. 
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Fig. 6. — Same as Figure 4 except that Mgrav is replaced by MdMi a^nd Mstars is calculated for Mstars/iy = 6.1 Mq / Lq. 



transition between stellar-dominated and DM-dominatcd 
regimes and yields Mstars/ Ly = Q-l Mq/ Lq consistent 
with the values determined in §6. 

In Figure 6 we plot Mgtars and Mdm corresponding to 
this fiducial stellar model. The inner data point cannot be 
seen in the figure because the dark mass is negative there, 
Mdm = -4.0 ± 1.7 x lO^M©, though it has large error. 
Consequently, we exclude the inner radial bin from fits to 
Mdm when using the fiducial value of Mstars/ Ly ■ 



We plot the best-fitting power-law and NFW models to 
-A^DM in Figure 6; the fit parameters are listed in Table 3. 
The quality of the fits and the parameter values of the DM 
models are consistent with those obtained for the gravitat- 
ing mass within the Itr errors. (Note: including the inner 
(negative) data point results in fits that are of significantly 
lower quality: — 19.4 for power law and = 18-5 for 
NFW, each with 3 dof.) Similar to the fits to Mgrav, the 
Moore model does not fit quite as well as NFW (x^ =3.6 
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for 2 dof), but the Hernquist model provides a comparable 
fit (x^ = 1.8 for 2 dof). The core parameter for the Hern- 
quist model is poorly constrained (a = 0.34 ± 0.77 Mpc). 

8. SYSTEMATIC ERROR BUDGET FOR MuM 

This section provides an assessment of the magnitude of 
systematic errors on the fits to Mdm- 

X-ray Background: To examine the sensitivity of the mea- 
sured dark matter properties to reasonable background er- 
rors we repeated our entire analysis using the standard 
background templates and correction files (§2) renormal- 
ized to have count rates ±20% of their nominal values. 
We find that the derived gas density profile is hardly af- 
fected while the temperature profile is modified within the 
IfT— 1.5(7 errors: p = (—0.09, 0.07) for these different back- 
ground levels. The effect on the DM fits is less than the 
1(7 errors in the DM model parameters as shown in Table 
4. We also examined whether the DM fits were system- 
atically different if the normalization of the background 
was allowed to be a free parameter during the fits. We 
found that the DM fits are fully consistent with the results 
quoted above. We note, however, that in this case the best- 
fitting temperatures 20%) and abundances 40%) 
are lower than the values listed in Table 1, although the 
differences are significant only at the « 2(7 level. 

Stellar Mass: In Table 4 we quote the ranges in model pa- 
rameters fitted to Mdm considering the variation in Mgtars 
implied by Mstars/^^y = (1 - 7.4) Mq/Lq (see §6). Note 
that for Msia.is/Lv = ^Mq/Lq the inner radial bin for 
Mdm is positive, and so we included it in the fits in Ta- 
ble 4. Uncertainties associated with the stellar mass are 
within the la errors. 

Center: Since the X-ray centroid varies with radius we ex- 
amined the sensitivity of our results to the centroid. For 
comparison to our results obtained using the centroid com- 
puted within a circular aperture of radius 20", we used a 
radius of 3'. Using this larger radius we obtain a centroid 
offset by « 8" from the previous value. However, we find 
that the results obtained for the DM fits change negligibly 
(< 1%) when using the different centroid. 

Deprojection Method: Our method to deproject the gas 
density and temperature uses the projections of param- 
eterized models (§4.2). For comparison, we also per- 
formed a non-parametric deprojection using the "onion- 
peeling" method as implemented in our previous studies 
(e.g., Buote 2000; Buote et al. 2003a: Lewis et al. 2003). To 
obtain stable deprojected parameters for the A2589 Chan- 
dra data we imposed the following restrictions: (1) the 
deprojected temperatures in annuli #3-5 were confined to 
the 1(7 limits in Table 1, and (2) the deprojected metallic- 
ities were set to those in Table 1. Because of these restric- 
tions the results obtained from this deprojection should 
only be viewed as an indicator of the sensitivity of the re- 
sults to deprojection and not as the preferred values. We 
found that the deprojected temperature and density pro- 
files give results consistent with those obtained from the 
parametric deprojection when fitted with the same mod- 
els in three dimensions. However, the gas density is more 
centrally peaked than indicated by the f3 model and thus, 
like in A2029 (Lewis et al. 2003), we find that adding a 
central cusp to the f3 model provides a better fit. We find 



that Mdm is small, but positive, in the central bin. Using 
all five data points we obtain the results reflected in Table 
4 for fits to the DM. 



9. CONCLUSIONS 

Our analysis of the Chandra data indicates that from the 
largest radius probed (r = 0.07rvir) down to r ~ 0.02rvir 
the dominant contributor to the gravitating matter in 
A2589 is the DM. Over this region the radial profile of 
the DM is fitted well by the NFW and Hernquist profiles 
predicted by CDM. (These models fit the gravitating mat- 
ter well down to r = 0.009rvir-) The inferred value of 
rvir of the NFW model is consistent (within the relatively 
large measurement errors) with that expected from CDM 
simulations (1.7 Mpc, see §4.3); equivalently, the concen- 
tration parameter c = 4.9 ± 2.4 is in the range expected 
for a massive cluster (Navarro et al. 1997). The DM pro- 
file over this region is also described well by a power law, 
M(< r) oc r'^ p (X r~" where a = 3 — 7. We find 
a = 1.35 ± 0.21 (statistical error, see §8) which is con- 
sistent with the NFW value a = 1 (and Moore value of 
a = 1.5) but is significantly larger than a w 0.5 found in 
LSB galaxies (e.g., Swaters et al. 2000) and expected from 
SIDM models (Spergel & Steinhardt 2000). We note that 
the profile of the gravitating matter (a = 1.37 ± 0.14) is 
also very consistent with the DM profile and with CDM 
predictions (El-Zant et al. 2003). 

The DM properties of A2589 agree very well with those 
obtained for A2029 (Lewis et al. 2003). Because the cores 
of these clusters are undisturbed by interactions from a 
central radio source, the results we have obtained for 
their core DM distributions provide critical confirmation 
of (in some respects) similar results obtained from other 
otherwise-relaxed clusters with disturbed cores (most no- 
tably Hydra-A, David et al. 2001) and clarify ambiguous 
cases like A1795 (Ettori et al. 2002). Furthermore, the 
evidence that the NFW-CDM profile is a good fit to clus- 
ter DM halos on larger scales (O.lrvir ^ < 0.5rvir) (e.g., 
Allen et al. 2002; Pratt & Arnaud 2002) is extended down 
to r ^ O.Olrvir by our complementary analyses of the 
Chandra data of A2029 and A2589. This agreement of 
cluster DM profiles with the NFW model is also supported 
by the recent weak-lensing analysis of several clusters by 
Dahle et al. (2003) and the joint lensing- X-ray study of 
Arabadjis et al. (2002). 

In sum, the current evidence from X-ray and lensing 
studies indicates that the radial DM profiles for cluster- 
sized DM halos are consistent with CDM predictions 
within the expected cosmological scatter (Bullock et al. 
2001). The disagreement between the NFW-CDM profile 
and the flatter proflles observed for LSB galaxies indicates 
that either CDM simulations are adequate for clusters but 
inadequate for LSB galaxies or the interpretations of the 
measurements of the DM profiles for LSB galaxies are in 
error. Some recent studies have indeed suggested that sys- 
tematic errors could account for the flat density profiles 
obtained for LSB galaxies (e.g., Swaters et al. 2003). 

Future X-ray and lensing studies of the DM in other (re- 
laxed) cluster cores are essential to refine the constraints 
on the DM properties obtained from these initial studies 
and to quantify precisely the cosmological scatter. 
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Table 4 

Error Budget for Dark Matter Fits 



Best AStatistical ABackground AMstars ADeproj 

Power Law: 

7 1.65 ±0.21 (+0.17,-0.16) (+0.06,-0.00) +0.35 
NFW: 

c 4.9 ±2.4 (+1.9,-1.7) (+1.2,-0.6) -3.8 

rvir(Mpc) 1.0 ±05 (+0.3,-0.2) (+0.06,-0.05) +1.9 



Note. — The "Best" column indicates the best-fit value and "AStatistical" the Icr statistical error from Tabic 3. "ABackground" gives 
the results when the X-ray background level is set to ±20% of nominal. "AMstars" represents the uncertainty associated with the stellar 
mass-to-light ratio. "ADeproj" provides an estimate of the error associated with the deprojection procedure. 
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